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Paradoxically enhanced cognitive processes in neurological disorders provide vital clues to understanding 
neural function. However, what determines whether the neurological damage is impairing or enhancing is 
unclear. Here we use the performance of patients with two disorders of the striatum to dissociate mech- 
anisms underlying cognitive enhancement and impairment resulting from damage to the same system. In 
a two-choice decision task, Huntington's disease patients were faster and less error prone than controls, 
yet a patient with the rare condition of benign hereditary chorea (BHC) was both slower and more error 
prone. EEG recordings confirmed significant differences in neural processing between the groups. Analysis 
of a computational model revealed that the common loss of connectivity between striatal neurons in BHC 
and Huntington's disease impairs response selection, but the increased sensitivity of NMDA receptors in 
Huntington's disease potentially enhances response selection. Crucially the model shows that there is a crit- 
ical threshold for increased sensitivity: below that threshold, impaired response selection results. Our data 
and model thus predict that specific striatal malfunctions can contribute to either impaired or enhanced 
selection, and provide clues to solving the paradox of how Huntington's disease can lead to both impaired 
and enhanced cognitive processes. 

© 2014 The Authors. Published by Elsevier Inc. 
This is an open access article under the CC BY-NC-ND license 
(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Paradoxically enhanced cognitive processes in neurological dis- 
orders both present stern challenges to our understanding of their 
underlying mechanisms and provide vital clues to furthering that un- 
derstanding (Frank et al., 2004; Shiner et al., 2012; Kapur et al., 2013). 
A key challenge is to solve the underlying paradox of how damage 
to the same neural system can result in both enhanced and impaired 
cognitive processes. We took a novel approach to tackling this chal- 
lenge by using the same cognitive task in two neurological disorders 
affecting the same system to dissociate their separate functional im- 
plications; we then used computational analyses to test the inferred 
hypothetical effects on the neural system. 



Abbreviations: AMPA, a-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid; 
BHC, benign hereditary chorea; EEG, electroencephalography; ERP, event related po- 
tential; GABA, -y-aminobutyric acid; MMN, mismatch negativity; NMDA, N-methyl- 
D-aspartate; RON, reorientation of attention; MSN, medium spiny neuron; FSIs, fast 
spiking interneurons; MMSE, Mini Mental Status Examination. 
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We previously showed that Huntington's disease patients have a 
paradoxical enhancement in a simple auditory decision-making task 
(Beste et al., 2008; Saft et al., 2008): on trials with irrelevant sensory 
information, a significantly enhanced mismatch negativity (MMN) 
signal from EEG recordings corresponded to significantly reduced er- 
rors and response times compared to controls. The MMN is defined 
as a phasic negativity of the event-related potential (ERP) of the EEG 
evoked by rare deviant stimuli in a sequence of stimuli and may 
reflect the recognition of rare events deviating from frequent back- 
ground events by the auditory system (Kujala et al., 2007; Naatanen 
et al., 2012). 

As Huntington's disease primarily affects the striatum, we hypoth- 
esized that such paradoxical enhancement of response selection is the 
result of the interaction between two factors in Huntington's disease 
progression (Beste et al., 2008): (i) increased NMDA receptor sensi- 
tivity at corticostriatal synapses (Fan and Raymond, 2007) and (ii) the 
consequent degeneration of medium spiny neurons through excito- 
toxicity, potentially driven by the increased calcium influx through 
the NMDA receptors. Yet, how or even if this interaction leads to an 
enhancement of striatal processing is unknown; and even if it does, 
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how this interaction then usually leads to impaired cognitive perfor- 
mance in Huntington's disease (Lawrence et al., 1998, 2000; Stout et 
al., 2001 ) is unknown. 

A critical test for these hypotheses would require an examination 
of sensory decision-making processes in a neurological condition that 
is characterized by medium spiny neuron dysfunction but without the 
corresponding enhancement of NMDA receptor function. In this re- 
gard, benign hereditary chorea (BHC, MIM 118700) may serve as a 
model (for review see Kleiner-Fisman and Lang, 2007 and Inzelberg 
et al., 2011). BHC is a rare autosomal dominant disease (prevalence 
between 1 and 2 in 1,000,000 people), which is caused by mutations 
in the TITF1 gene on chromosome 14ql3 encoding the thyroid tran- 
scription factor-1 (also known as TITF1, TEBP or NKX2-1) (Inzelberg 
et al., 201 1 ). TTF1 mediates striatal interneuron migration from the 
medial ganglionic eminence (a precursor of the globus pallidus) to 
the lateral ganglionic eminence (a precursor of the striatum) (Sussel 
et al., 1999). This migration is reduced in TITF1 mutation and leads 
to dysgenesis of the formation of medium spiny neurons (Kleiner- 
Fisman et al., 2005; Yoshida et al., 2012). By contrast, TITF1 has no 
known role in NMDA receptor function (Sussel et al., 1999) and thus 
NMDA receptor function in BHC is likely unaltered. 

Consequently, we tested a BHC patient and manifest Huntington's 
disease patients against age-matched controls to dissociate the con- 
tributions of striatal changes in NMDA sensitivity and network struc- 
ture to enhanced response selection. Using a computational model of 
the striatum, we tested the hypothesis that either of these changes in 
the striatum could directly affect response selection, and then tested 
the hypothesis that increased NMDA sensitivity in the striatum of 
Huntington's disease patients could enhance selection despite the 
loss of neurons. Critically, we sought from the model the key expla- 
nation for the paradox of why both impaired (Jahanshahi et al., 1993; 
Lawrence et al., 1998) and enhanced response selection could arise in 
Huntington's disease. 

2. Materials and methods 

2.1. Patients and participants 

A group of manifest Huntington's disease patients (N = 10; 6 
females) between 20 and 40 years of age (mean age 33.45 ± 5.5) 
without any medication and a group of controls (N = 10; 6 females) 
matched to the manifest Huntington's disease group in age, sex and 
educational background were recruited. Along with these groups, a 
case of benign hereditary chorea (BHC) (female, 24 years of age) with 
strong motor symptoms was recruited and compared to the other 
groups using single-case t-statistics (Crawford and Garthwaite, 2012). 
BHC was genetically confirmed by detecting the mutation in the TTF1 
gene and the BHC case did not take any medication at the time point 
of examination. The clinical and neuropsychological test data of these 
groups and the BHC case are shown in Table 1. The study was ap- 
proved by the Ethics Committee of the Ruhr-University of Bochum. 
The study was conducted according to the Declaration of Helsinki. All 
participants and the BHC case gave written informed consent. 

2.2. Task 

Subjects performed a distraction paradigm (Beste et al., 2008; see 
also Schroger and Wolff, 1998) in which tones at three different fre- 
quencies ( 1 000 Hz, 1 1 00 Hz, 900 Hz) were presented for either 400 or 
200 ms. One pitch (i.e., 1000 Hz) served as the standard tone, which 
was presented in 80% of trials. The other pitches were presented with 
a frequency of 10% each. The subjects were asked to respond with 
their thumb and indicated, whether the tone was short (right button 
press) or long (left button press). Variations in the pitch of the tone 
thus served as distraction. 



2.3. EEC recording and analysis 

The EEC was recorded from 32 scalp electrodes at standard po- 
sitions according to the 10/20 system. Electrode signals were sam- 
pled at 1000 Hz with Cz as primary reference. The resulting time- 
series were downsampled to 256 Hz in offline post-processing and 
re-referenced to linked mastoids. After a first visual inspection of 
the time-series a bandpass filter from 0.5 to 20 Hz (48 dB/oct) was 
applied. Ocular artefacts (blinks and saccades) as well as pulse arte- 
facts were corrected using independent component analysis (Info- 
max algorithm), applied to the un-epoched time-series. Besides the 
MMN, other processes related to the reorientation of attention (RON 
(Schroger et al., 2000)) have also been found to be increased in their 
efficacy in Huntington's disease, whilst attentional shifts (reflected by 
the P3a; Escera and Corral, 2007) are not affected (Beste et al., 2008). 
If exaggerated glutamatergic neural transmission is of similar im- 
portance for reorientation processes, the RON and P3a mechanisms 
should be similarly modulated to the MMN. To measure the MMN, 
P3a and the RON the EEC time-series were epoched in segments from 
-200 till 800 ms after the stimulus onset. Within these epochs an au- 
tomated artefact rejection procedure was applied. Rejection criteria 
included a maximum voltage step of more than 60 (iV/ms, a maxi- 
mal value difference of 150 u.V in a 250 ms interval or activity below 
0.1 |xV. After this, a baseline correction was -200 till 0 (i.e., time point 
of stimulus presentation was applied). To measure the MMN, P3a and 
RON, difference waves were calculated (distractor minus standard 
ERPs) (Kujala et al., 2007). In these difference waves the MMN was 
defined as the most negative peak between 100 and 250 ms. The P3a 
was defined as the most positive peak between 250 and 500 ms and 
the RON was defined as the most negative peak between 400 and 
600 ms post-stimulus presentation. The ERPs and the behavioral data 
are shown in Fig. 1 . 

2.4. Statistical tests 

Behavioural data of the BHC case in relation to the Hunting- 
ton's disease group and controls were analysed using single-case t- 
statistics (Crawford and Garthwaite, 2012) applying Crawford and 
Howell's method. This method is widely used and offers the best pos- 
sible way to compare single cases with groups of other subjects (for 
review see Crawford and Garthwaite, 2012). ANOVAs were used to ex- 
amine effects between the Huntington's disease and control groups. 
Post-hoc tests in the Huntington's disease and control group were 
Bonferroni-corrected where necessary. All descriptive statistics are 
given as mean and standard error of the mean. 

2.5. Population model of competing striatal populations 

We used reduced models of the striatal microcircuit (Figs. 2 and 3) 
to understand the mechanisms by which Huntington's disease- and 
BHC-like changes in the striatum could potentially alter selection. In 
the baseline model, each medium spiny neuron population's firing 
rate a was modelled as a leaky integrator of inputs from the cortex 
and other medium spiny neuron population: 

raj = -cii + Wji|_ajJ + + wilj (1) 
r a, = -q + wy LdjJ + + w^j (2) 

which decays with time constant T, and is driven by the sum 
of weighted cortical input Ij.Ij and feedback inhibition (Wy.Wjj < 0) 
from the other medium spiny neuron population. The operation |_-J + 
indicates that negative medium spiny neuron population activity is 
rectified to zero. 



C. Beste et al. / Neurolmage: Clinical 4(2014) 623-634 



625 




Fig. 1. Behavioural and electrophysiological differences in the auditory discrimination task. (A) Error rates and reaction times on distractor (black) and standard trials (white) in 
controls, manifest Huntington's disease and the BHC case ( ± SEM). The data shows better selection performance (fewer errors, faster reaction times) in Huntington's disease, but 
reduced selection performance in BHC, compared to controls. (B) Differences (standard minus distractor) at electrode Fz for controls (red), manifest Huntington's diseases (black) 
and BHC (blue). The different parts of the difference waves are labelled MMN, P3a and RON. Paralleling the behavioural data, the MMN and RON were increased in Huntington's 
disease and decreased in BHC, relative to controls. 

Table 1 

Description and analysis of performance (test scores) in different neuropsychological tests. 



Test BHC HD Controls HD vs. controls BHC vs. controls BHC vs. HD 



UHDRS motor 




26 


24.5 (9.9) 


NA 






score 














CAC repeat 




NA 


45.3 (5.5) 


NA 






length 














IQ(MWT-B) 




105 


103 (4.5) 


113(8.7) 






Stroop test 


Colour naming 


61 


48.5(16.7) 


81.1 (3.6) 






(UHDRS) 
















Colour reading 


9S 


65.6(14.9) 


93.7(5.1) 








Interference 


36 


30.5(12.5) 


52.1 (3.5) 






Symbol-digit test 




51 


39.4(10) 


55.3 (3.4) 






(WA1S) 














Word fluency 




39 


26(8.8) 


46.2 (9.6) 






(Benton) 














Digit span 


Forward 


8 


5.1(1.1) 


9.1(1.1) 




** 


(WMS-R) 
















Backward 


5 


2.2 (0.9) 


8.0(0.9) 






Block span 


Forward 


8 


4(0.7) 


8.7 (0.9) 






(WMS-R) 
















Backward 


7 


3.1 (0.5) 


7.9 (0.9) 






Benton test 




11 


7.6(1.8) 


13.5(0.9) 






(visual memory) 














Mini Mental 




30 


28.5(2.1) 


30 


*** 


** 


Status 














Examination 














(MMSE) 















Significances denote: "p < .01. *"p < .001. 



Extending this model to include the self-inhibition of each medium 
spiny neuron population gave: 

ra, = -aj + WjiLajJ + + wilj + w,iLaiJ+ (3) 
r a, = -Oj + Wjj |_aij + + wilj + Wjj |_QjJ + (4) 

where Wjj.Wjj are the self-inhibition weights. This model is 
schematically illustrated in Fig. 3B. 

Further extending this model to incorporate the feedforward in- 
hibition of the medium spiny neuron populations by a population of 
fast spiking interneurons gave: 

rdj = -dj + WjiLajJ + +wilj + w,iLaiJ+ + w F I F (5) 

TOj = -Oj + Wjj LdjJ + + W[lj + Wjj |_OjJ + + w F l F (6) 



where 1 F = f(Ii,Ij) is the input received from the fast spiking in- 
terneuron population, and is some (strictly-increasing) function of 
the cortical inputs to the striatum. Given the dense connectivity of 
each fast spiking interneuron with their target medium spiny neu- 
rons (Humphries et al., 2010; Planert et al., 2010), we have assumed 
here a "broadcast" model in which the fast spiking interneuron pop- 
ulation projects equally to the two competing medium spiny neuron 
populations with weight w F . This model is schematically illustrated in 
Fig. 3C. We note that fast spiking interneurons do not express NMDA 
receptors and so are likely spared in Huntington's disease. 

Inputs from the cortex emulated a pre-selection and a selection 
phase of a two-choice task; the success of selection was assessed 
in the selection phase. In the pre-selection phase both response- 
representing medium spiny neuron populations received from the 
cortex the same input r pre , simulating the onset of task-salient infor- 
mation. In the selection phase, the first medium spiny neuron pop- 
ulation received a step in input of r ste p, simulating an increase in 
salience of that response after the perceptual decision was made (e.g. 
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the pushing of the left or right button after the decision on whether 
the tone was short or long). 

The consequent output of the two medium spiny neuron popu- 
lations was read-out in a signal selection framework (Curney et al., 
2001): we considered that the response was correctly selected if the 
output of the "winning" response population crossed a threshold 0 H 
above its steady-state and that the output of the "losing" response 
population crossed a threshold 6i below it steady-state; otherwise 
we consider the selection to be ambiguous (the formal expression 
is given by Equations 6-7 of the Supplementary results). Here we 
assume that ambiguous response selection leads to both a potential 
error and an increased reaction time as the ambiguity must be instead 
resolved in the final common motor pathway to produce a motor out- 
put (Prescott et al., 1999). For all illustrative plots (Figs. 2-5) we use 
upper and lower thresholds of ± 2 Hz. 

The D1-D2 model (Fig. 5) splits the two medium spiny neu- 
ron populations into two sub-populations of Dl and D2-expressing 
medium spiny neurons. The dynamical equations of each popula- 
tion were identical to those for the two-population models, except 
for the addition of input from all three other medium spiny neu- 
ron sub-populations (schematically illustrated in Fig. 5). This four- 
dimensional system did not allow for exact solution, and so results 
were obtained by numerical integration using the Euler numerical 
method. The time constant was arbitrarily set to t = 1 , as it has no ef- 
fect on the equilibrium state. To enable direct comparison, all weights 
were set to the same values as those used to generate the illustrative 
analytical results in Fig. 3: lateral and self-inhibition weight = -0.5; 
fast spiking interneuron input weight = -0.1. 

2.6. Biophysical model of competing striatal populations 

We examined whether or not the introduction of non-linear neu- 
ronal dynamics would alter the conclusions from the reduced model 
by constructing a canonical microcircuit model. Our biophysical 
model replicated the reduced fast spiking interneuron feed-forward 
inhibition model as this was the most complete model we examined. 
The biophysical model also provided a strong test of the striatum's 
ability to perform signal selection as it explicitly incorporated the 
weak synaptic efficacy of medium spiny neuron local collateral con- 
nections and the realistic connection densities of the striatal circuit. 

Two model medium spiny neurons represented the output of two 
small medium spiny neuron populations encoding two alternative 
responses. They received spike train input from two cortical popu- 
lations encoding the salience of the alternative responses. The other 
components of the model simulated spiking input from the other 
medium spiny neurons and from the fast spiking interneurons of the 
striatum (Tepper et al., 2004). 

We used the spiking medium spiny neuron model from Humphries 
et al. (2009a), which reproduces characteristic behaviour of these 
neurons including their f-1 curves, slow-rise to first spike on current 
injection, and paired-pulse facilitation. AMPA and GABAa induced 
post-synaptic currents were modelled as standard single exponential 
functions of short decay. As NMDA receptors are primarily affected 
in Huntington's disease, we used a detailed NMDA synaptic current 
model including voltage-dependent block by magnesium ions and 
saturation (full equations and parameters of the biophysical model 
neuron are given in the Supplementary methods). Enhanced NMDA 
receptor sensitivity was modelled by increasing the NMDA receptor 
conductance. 

Each of the cortical, fast spiking interneuron, and medium spiny 
neuron populations were simulated as a generator of spiking input to 
the model medium spiny neurons, producing the spike-events that 
occur across N afferents to the medium spiny neuron. Given time-step 
Af (s) and mean spike rate r (spikes/s), the probability of a spike per 
afferent is p(s) = rAt. The total number of spike-events S at each time- 
step is then just drawn from a binomial distribution S = B(N,p(s)). The 



resulting time-series of spike-events is equivalent to the pooling of N 
spike trains modelled as independent renewal processes. 

For the cortex, two generators were used with mean rates rj and 
r 2 . Following arguments in Humphries et al. (2009b) for the size of 
the active afferent cortical population to one medium spiny neuron 
we set N = 500 for both. For the feed-forward fast spiking interneuron 
population, the mean rate was based on the rate of cortical firing. We 
simulated our large-scale network model of the striatum (Humphries 
et al., 2010) with a range of cortical inputs to determine the fast 
spiking interneuron population's input-output function. When sim- 
ulating the biophysical microcircuit model, we took the mean of the 
cortical inputs and l^and looked up the corresponding mean fast spik- 
ing interneuron rate from the derived input-output function. We 
assumed a density of 1% fast spiking interneurons in the striatum, 
which gives N = 31 fast spiking interneurons afferent to 1 medium 
spiny neuron (Humphries et al., 2010). 

The feedback medium spiny neuron population's mean rate was 
based on the output of the two neuron models. We took the mean 
of their spike-rate over a 1 s window. The medium spiny neuron 
feedback was turned on after 1.5 s to allow stable firing to develop. 
We used N = 728 medium spiny neurons afferent to 1 medium spiny 
neuron (Humphries et al., 2010). Excitotoxic death of medium spiny 
neurons in Huntington's disease and loss of medium spiny neuron 
local fibres in BHC were modelled as a reduction in the size of this 
feedback population. 

The task-emulating inputs to the biophysical model had the same 
protocol as those for the reduced model (see above), with the only 
difference being that the cortical input was now explicitly given as 
spiking input. Here we used r pre = 2 Hz (giving 1000 spikes/s total 
input to each medium spiny neuron). The step was made at 7 s, and 
firing rate estimates were made in 5 s windows before and after the 
step. The long durations in the biophysical model were used only 
to allow accurate assessment of firing rates from the resultant spike 
trains of the two model neurons. 

As reported in the Results section, we searched for the minimum 
step-size that would achieve unambiguous signal selection. For each 
increment of the search, we increased the input step by 0.05 spikes / s 
per cortical train, and repeated the full simulation 20 times. Across 
these simulations we computed the mean firing rate in the post-step 
window for both the winning and the losing medium spiny neuron 
population, and compared to the upper and lower thresholds. The 
search was halted once both thresholds were crossed. When com- 
paring minimum step-sizes between models (Fig. 4B), we considered 
them different if there was a greater than ± 1 SEM overlap between 
the mean values. 

3. Results 

3. 1. Striatal disorders have opposite effects on task performance 

We used a simple auditory decision task to probe the effects of 
striatal disorders on response selection. On every trial subjects were 
asked to indicate which of a short (200 ms) or long (400 ms) tone was 
played; 20% of trials had an irrelevant change in pitch as a distractor. 
The task was carried out by a group with manifest, non-medicated 
Huntington's disease (n = 10), a matched control group (n = 10), and 
a single BHC patient. Compared to the control group and BHC patient, 
the Huntington's disease group was impaired on a range of standard 
cognitive tasks (Table 1 ). 

The Huntington's disease and control groups were newly tested, 
and so also acted as a stringent test that we could replicate our 
previous surprising results showing enhanced response selection in 
Huntington's disease patients (Beste et al., 2008). We found that 
the three groups differed in their error rates, but only on distrac- 
tor trials (Fig. 1). Error rates were generally higher for distractor 
trials (7.4 ± 0.47%) compared to standard trials (3.93 ± 0.11%) 
(F(l,18) = 41.00; p < 0.001; =0.695). 
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Fig. 2. Signal selection in healthy-state, Huntington's disease-like and BHC-like models 
of the striatum. The medium spiny neuron (MSN ) populations receive inputs represent- 
ing competing response signals, with the more salient response having a step increase 
(top); response selection is read-out from the population outputs after the step, by 
checking that both upper (0 H ) and lower (8i) thresholds are crossed. All example sim- 
ulations show responses to the same input. (A) The healthy-state model. Line-widths 
in the model schematic (right) indicate relative synaptic weights. (B) A BHC-like model, 
with loss of medium spiny neuron intra-connectivity. (C) A Huntington's disease-like 
model, with loss of medium spiny neuron intra-connectivity due to cell death and 
enhanced cortical input weight due to increased NMDA receptor sensitivity. 



The mixed effect ANOVA (standard vs. distractor as within-subject 
factor, Huntington's disease vs. control as between subject factor) re- 
vealed an interaction "standard/distractor x group" (F( 1,1 8) = 14.19; 
p = 0.001 ; =0.441 ). Post-hoc tests showed that the control and Hunt- 
ington's disease groups did not differ on standard trials (p > 0.6) 
but did on distractor trials (p < 0.001), with Huntington's disease 
patients showing better performance compared to controls (cf. Fig. 
1 A). Importantly, the BHC patient was also less accurate on distractor 
trials, when compared to the control group (t = 2.4; p = 0.019) and 
Huntington's disease group ( t = 2.6; p < 0.01 ). No difference between 
the BHC patient and the control or Huntington's disease groups was 
evident in standard trials (p > 0.4). 

We found that the three groups also differed in their reaction 
times. A mixed effects ANOVA (standard vs. distractor as within- 
subject factor, Huntington's disease vs. control as between subject 
factor) showed that reaction times were faster on standard trials 
(574 ms ± 10.1) compared to distractor trials (613 ± 11.1 ms) 




-0.4 -0.2 
Inhibition weight 



Fig. 3. Selection by the striatum is differentially modulated by increased NMDA sensi- 
tivity and loss of striatal connectivity. In each panel, we show a schematic of the model 
circuit (right), and the analytically-computed minimum input step required for unam- 
biguous signal selection (left) as a function of the inhibitory and input weights. Each 
heat map's colour intensity encodes the required minimum input step (firing rate, in 
Hz) for those weights, and the white line is the analytical boundary that separates the 
regions of smaller (above) and larger (below) step-sizes than the healthy-state model 
(bottom left-hand corner). Heat maps are plotted for specific choices of model parame- 
ter values to illustrate the general insights from the analytical solutions (Supplementary 
results), and so these results do not depend on specific choices of parameter values. (A) 
The baseline model of weakly-competing medium spiny neuron (MSN) populations. 
(B) The effect of adding inhibitory feedback to competing medium spiny neuron popu- 
lations. (C) The effect of adding feed-forward inhibition from fast spiking interneurons 
(FSIs) to competing medium spiny neuron populations. The fast spiking interneuron 
population is a "broadcast" model that equally samples from the cortical inputs and 
equally outputs to the medium spiny neuron populations. 

A 



B 




Missing MSNs (%) 

Fig. 4. Tradeoff between increased NMDA sensitivity and decreased striatal connec- 
tivity is robust to nonlinear neuronal dynamics. (A) The reduced spiking microcircuit 
model. Two spiking model medium spiny neurons (red, blue) represent the two pop- 
ulations of medium spiny neurons (MSNs). Each neuron model receives glutamatergic 
input (arrows) via AMPA and NMDA receptors from its afferent cortical populations, 
and GABAergic input (circles) from striatal fast spiking interneurons (FSIs) and from 
the other medium spiny neurons. Each box represents a source of spikes based on the 
number of spikes received from the indicated sources. The numbers of connections 
between each indicated population are derived from a full three-dimensional model of 
striatal micro-anatomy (Humphries etal., 2010)). Inset: an example of the output of the 
two medium spiny neuron models during a single selection test. Cortical input to MSN1 
increased at 7 s (grey line), representing a more salient response signal. (B) Landscape 
of minimum input step required for unambiguous signal selection (in total spikes/s 
from cortical spike-trains). The value for each parameter pair (missing medium spiny 
neurons, NMDA conductance) is the mean of 20 numerical searches for the minimum 
input step. The white line gives the boundary between better (smaller step) and worse 
(larger step) signal selection performance than the healthy-state model (bottom left). 
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Fig. 5. Tradeoff between increased NMDA sensitivity and decreased striatal connectivity occurs for selective loss of D2 medium spiny neurons (MSNs). (A) The Dl -D2 microcircuit 
model. Right: we split the two medium spiny neuron populations into separate sub-populations of medium spiny neurons predominantly expressing the Dl and D2 receptors 
for dopamine, and modelled the resulting connections within and between them. Both the Dl and D2 medium spiny neuron populations received the same input signal from 
their corresponding cortical population; the selection was read-out from the Dl medium spiny neuron population's output. Left: minimum input needed to achieve unambiguous 
signal selection following selective loss of the D2 medium spiny neurons only: all D2 medium spiny neuron output weights were changed whilst the Dl medium spiny neuron 
output weights were held constant. Heat map conventions are the same as in Fig. 3; the results plotted here are from numerical simulations. (B) The D1-D2 microcircuit model 
with the addition of feedback self-inhibition to each medium spiny neuron population: the heat map gives the minimum input needed to achieve unambiguous signal selection 
following selective loss of the D2 medium spiny neurons only. (C) The D1-D2 microcircuit model with the further addition of feed-forward fast spiking interneuron (FSI) input to 
each medium spiny neuron population: the heat map gives the minimum input needed to achieve unambiguous signal selection following selective loss of the D2 medium spiny 
neurons only. 



(F(l,18) = 7.62; p = 0.013; =0.297). There was no interaction 
"standard/distractor x group" (F(l,18) = 0.03; p > 0.9), but the 
main effect group revealed that the Huntington's disease group 
showed faster responses (550 ± 11.8 ms) than the control group 
(637 ± 12.1 ms) (F(l,18) = 27.02; p < 0.001; =0.600). For the BHC 
patient, the single-subject t-tests (see Materials and methods) re- 
vealed that in distractor trials reaction times were longer compared 
to those in both the control group (t = 3.1 1 ; p < 0.006) and the Hunt- 
ington's disease group (t = 5.26; p < 0.001 ), but not in standard trials 
(p > 0.5). 



3.2. Performance changes have specific neural correlates in EEG 

The electrophysiological data (difference wave ERPs) (MMN, P3a 
and RON) for the control group, Huntington's disease group and BHC 
patient are shown in Fig. IB. The mismatch negativity (MMN) was 
analysed at electrode Fz (e.g. Naatanen et al., 2012). Here, the Hunting- 
ton's disease group revealed a larger MMN (-9.99 ± 0.66 |xV) com- 
pared to the control group (-7.32 ± 0.45 ii.V) (univariate ANOVA; 
F(l,18) = 17.68; p = 0.001; =0.496). The latency of the MMN was 
shorter in the Huntington's disease group (119 ± 6.5 ms) compared to 
the control group (152 ± 5.5ms)(F(l,18)=14.63;p=0.001;=0.448). 



Importantly, the BHC patient revealed a smaller MMN (-3.12 jiV) 
compared to the Huntington's disease group (f = 4.51; p < 0.001) 
and the control group (t = 2.81; p = 0.01). In the BHC patient, the 
MMN latency (160 ms) was longer compared to that in the Hunt- 
ington's disease group (119 ± 6.5 ms) (t = 2.17; p = 0.02), but not 
different to controls (p > 0.4). For the P3a measured at electrode 
Fz, there was no difference between the Huntington's disease group, 
control group and BHC patient (all p > 0.5). For the RON, however, 
the pattern of results was identical to that in the MMN data. For 
the RON, the amplitude was larger (most negative) for the Hunting- 
ton's disease group (-7.85 ± 0.6 compared to the control group 
(-4.82 ± 0.61 iiV) (F(l,18) = 12.83; p = 0.002; =0.416). The RON 
in the BHC patient (0.2 ii.V) was less negative compared to that in 
both the control group (t = 2.86; p = 0.009) and Huntington's disease 
group (t= 3.67; p = 0.002). 



3.3. Hypotheses for response selection in the striatal microcircuit 

The disorders caused opposite changes in response selection com- 
pared to controls: reduced errors and response time for Huntington's 
disease patients (replicating our previous study, Beste et al., 2008); 
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increased errors and response time for the BHC patient. The electro- 
physiological data showed significant differences in event related po- 
tentials between the Huntington's disease group, control group, and 
BHC patient, confirming that the performance differences reflected 
changes in neural processing. 

We propose that this dissociation between Huntington's disease 
and BHC effects on response selection supports two hypotheses. First, 
that as both disorders predominantly affect the striatum, so com- 
putational processes involving the striatum directly influence the ob- 
served response selection; it thus follows that changes in the striatum 
are the cause of changes in this selection. Second, that as both disor- 
ders cause changes to the striatal network, the selection enhancement 
seen in Huntington's disease patients is solely due to enhanced NMDA 
sensitivity in the striatum. 

We simulated a striatal microcircuit model to test these hypothe- 
ses. To address the first hypothesis we sought to show how selec- 
tion could be either impaired or enhanced at the level of the striatal 
microcircuit itself. We further sought from the model a potential ex- 
planation for the key challenge: why, if our hypotheses are correct, 
Huntington's disease patients are usually impaired compared to con- 
trols in other decision-making tasks. Consequently, we focussed here 
on modelling the generic striatal mechanisms for response selection 
potentially common to a range of tasks. 

3.4. Computational analysis of enhanced and impaired response 
selection 

We first show analytically how Huntington's disease-like and BHC- 
like changes affect response selection in a reduced model of weakly- 
competing medium spiny neuron populations. We then show the 
robustness of these results to variations in the striatal circuitry, to the 
addition of biophysical details that may affect selection processes, and 
to the selective loss of striato-pallidal cells in Huntington's disease. 

3.5. Enhanced NMDA sensitivity compensates for depletion of striatal 
connectivity 

The computational role of the striatum is contentious. The pres- 
ence of local collateral connections between the GABAergic medium 
spiny neurons seems consistent with a winner-take-all computation 
on its inputs (Beiser et al., 1997), but the winner-take-all hypothesis 
has been undermined by the weakness of the synapses formed by 
these local collaterals (Jaeger et al., 1994; Czubayko and Plenz, 2002; 
Tunstall et al., 2002; Taverna et al., 2004). However, winner-takes-all 
is an extreme of a continuum of selection computations (Fukai and 
Tanaka, 1997). A more general form is signal selection, where an un- 
ambiguous selection between competing signals is indicated by the 
winning signal rising above an upper threshold and the losing sig- 
nals) falling below a lower threshold. The existence of many weak 
medium spiny neuron inputs to each medium spiny neuron is poten- 
tially sufficient to collectively influence its output (Humphries et al., 
2010; Chuhma et al., 201 1 ; Yim et al., 201 1 ). We also note that, as the 
primate striatum is predominantly silent (Hikosaka et al., 1 989), to see 
competition between medium spiny neuron populations would need 
a task invoking strong cortical inputs to those populations. Under 
these conditions, even weakly competitive medium spiny neurons 
are capable of modulating signal selection (Supplementary results, 
Section 1.1). We show here that this signal selection is modulated by 
Huntington's disease-like and BHC-like changes. 

We constructed a reduced model of two weakly competitive 
medium spiny neuron populations representing alternative responses 
(Materials and methods). Input to the model represents the putative 
cortical signals for the salience of the two potential responses (short/ 
long) on one trial: at the onset of the auditory tone, both inputs are 
set to some non-zero value, indicating that only these responses are 



salient. A further step in input to one medium spiny neuron pop- 
ulation then represents the decision-point at which the short/long 
percept is resolved (Fig. 2, top), and indicates the increased salience 
of the corresponding response. 

The output of the two medium spiny neuron populations after the 
step in input was read-out for signal selection (Fig. 2). Fig. 2A shows an 
example of signal selection by the medium spiny neuron population 
output for the healthy-state model. We compared the performance 
of this model to BHC-like and Huntington's disease-like models that 
emulated the changes to the striatal microcircuit thought to occur in 
these disorders. The loss of inter-medium spiny neuron connectivity 
in Huntington's disease (due to medium spiny neuron cell death) and 
BHC (due to depletion of medium spiny neuron local collaterals) was 
emulated by reducing the lateral inhibition weights between the two 
populations (Fig. 2B). The increased NMDA sensitivity in Huntington's 
disease was emulated by increasing the weight of input from the 
cortex to the two medium spiny neuron populations (Fig. 2C). Fig. 2B 
and C illustrates how example BHC- and Huntington's disease-like 
changes altered the signal selection performance compared to the 
healthy-state model. 

To understand the impact of these Huntington's disease- and BHC- 
like changes, we quantified response selection performance by finding 
the minimum step change (Al) in cortical input needed to get unam- 
biguous signal selection. This value thus indicates the difficulty of 
the selection process, and we could then quantify how that difficulty 
changes with alterations to the microcircuit. We exactly solved the 
models to find the minimum input step required in the healthy-state 
model, and after all possible BHC-like and Huntington's disease-like 
changes to the microcircuit (Supplementary results, Section 1 ). 

We find that, given the two thresholds (6> H . #l) for unambiguous 
signal selection, the minimum step in input to reach both is 

zj i = e L (i-w ijWji ) ( 7) 

WlWjj 

where wi is the cortical input weight, and (wy.Wjj) are the lateral 
inhibition weights between the two medium spiny neuron popula- 
tions (note that only the lower threshold # L appears in this function, 
as the step-size needed to reach that threshold is always at least as 
big as that needed to reach the upper threshold — see Supplementary 
results, Section 1.2). 

Eq. (7) shows that loss of connectivity between medium spiny 
neuron populations, reducing the magnitude of Wjj and Wjj, always 
increased the minimum step size, indicating impaired selection com- 
pared to the healthy-state model. This is consistent with the BHC pa- 
tient's poor selection performance compared to controls. Conversely, 
Eq. (7) shows that solely increasing the cortical input weight (wi) de- 
creased the minimum step size, indicating comparatively enhanced 
selection compared to the healthy-state model. The heat map in Fig. 
3A illustrates these main analytical insights by plotting the depen- 
dence of the minimum input step (in Hz) needed for unambiguous 
signal selection on the reduction in lateral inhibition weight and the 
increase in cortical input weight; the healthy-state model is given at 
the lower-left corner. 

In the Huntington's disease-like case, both decreased medium 
spiny neuron connectivity and increased cortical input weight occur 
simultaneously. Given some change in medium spiny neuron popula- 
tion connectivity from (wy.Wjj) in the healthy-state to (w?j, wjj) in the 
Huntington's disease-like state, we can determine from the model 
(Supplementary results, Section 1.2) the value of the cortical input 
weight wf needed to exactly compensate and maintain the same se- 
lection performance: 

WlW.j (l - w*-wt) 
wf= ^ (8) 

W M 0 - W :, W ) 



630 



C. Beste et al. / Neurolmage: Clinical 4(2014) 623-634 



We plot an example of this function as the white line in the heat 
map of Fig. 3A. We can immediately see that Eq. (8) defines a criti- 
cal threshold for increased NMDA sensitivity that separates states of 
"enhanced" and "impaired" selection compared to the healthy-state 
model. In the enhanced state, sufficiently increased NMDA sensitivity 
can still result in enhanced selection despite loss of medium spiny 
neuron inter-connectivity (region above the white line in the heat 
map of Fig. 3A). 

This is consistent with the Huntington's disease patients' en- 
hanced selection performance compared to controls. The analytical 
results (Eqs. (7) and (8)) show that the changes in selection due to 
changes in lateral inhibition and input weights, and therefore also 
this trade-off between them, exist irrespective of particular choices 
of threshold or weights. Therefore, our reduced model provides evi- 
dence for both our hypotheses: first, that the striatum can implement 
a selection process and thus directly influence response selection; 
second, that selection enhancement seen in Huntington's disease pa- 
tients is solely due to enhanced NMDA sensitivity. 

Moreover, our analytical results show that the impaired state also 
must always exist, in which insufficiently increased NMDA sensitiv- 
ity will result in impaired selection through the loss of medium spiny 
neuron intra-connectivity (region below the white line in the heat 
map of Fig. 3A). The existence of this region is crucial: after all, in most 
simple decision-making tasks Huntington's disease patients' perfor- 
mance is significantly worse than control subjects (Lawrence et al., 
1998), and so a model predicting a uniformly improved performance 
through increased NMDA sensitivity would be inconsistent with the 
proposed role of the striatum in response selection (implications of 
this are considered further in the Discussion section). 



3.6. Trade-off of NMDA sensitivity and connectivity loss is robust to 
details of the striatal circuit 

Our model of weakly-competing medium spiny neuron popula- 
tions provided evidence for our two hypotheses, but other sources 
of inhibition of medium spiny neurons within the striatum raised 
two open questions: first, whether they prevented signal selection 
from occurring at all; second, even if signal selection was maintained, 
whether the Huntington's disease- and BHC-like models still showed 
signal selection changes that were consistent with the behavioural 
data. 

We first examined the role of self-inhibition by each medium spiny 
neuron population. Local collateral connections in the above model 
are restricted to those between response-representing populations. 
Whilst this provides a baseline model for quantifying the effect of 
mutual competition on response selection, there is no a priori reason 
to suppose that medium spiny neurons in the two populations selec- 
tively connect only to medium spiny neurons in the other population: 
thus it is likely that medium spiny neurons within each population 
are also connected by local collaterals and consequently each popu- 
lation is self-inhibiting. We thus implemented a model incorporating 
feedback inhibition of each medium spiny neuron population (Fig. 
3B). 

We found that adding such self-inhibition does not alter the qual- 
itative effects of losing connectivity between medium spiny neu- 
rons and of enhanced NMDA sensitivity (Supplementary results, Sec- 
tion 2.1). As an illustration of these results, Fig. 3B shows that self- 
inhibition increases the minimum input step needed compared to 
the baseline model, but the critical threshold separating model states 
with enhanced and impaired selection remains (Supplementary re- 
sults, Eq. 16). Consequently, even with self-inhibition of the response- 
representing populations, BHC-like changes always lead to impaired 
selection, and Huntington's disease-like changes can lead to enhanced 
selection, again providing evidence for our two hypotheses. 



Short time-scale dynamics in the striatum are strongly depen- 
dent on feed-forward input from GABAergic fast-spiking interneu- 
rons to the medium spiny neurons (Tepper et al., 2004; Mallet et 
al., 2005;Humphries et al., 2009b, 2010). Thus we also examined its 
potential influence on the striatal circuit's response selection ability, 
and how that ability is altered by Huntington's disease- and BHC-like 
changes. Given the extensive collateralization of the fast spiking in- 
terneuron axon and the vast number of medium spiny neurons each 
fast spiking interneuron contacts (Humphries et al., 2010; Planert et 
al., 2010), there is no a priori reason that a population of fast spik- 
ing interneurons would preferentially connect to one local medium 
spiny neuron population . Thus we implemented a "broadcast" model 
in which a population of fast spiking interneurons projects in equal 
proportion to both the medium spiny neuron populations represent- 
ing alternative responses (Fig. 3C). 

We found that introducing such feed-forward inhibitory input 
does not alter the qualitative effects of losing connectivity be- 
tween medium spiny neurons and of increased NMDA sensitivity 
(Supplementary results, Section 2.2). As an illustration of these re- 
sults, Fig. 3C shows that feed-forward inhibition from the fast spiking 
interneurons decreases the minimum input step needed compared to 
the baseline model (see Supplementary results, Section 2.2 for further 
discussion), but the critical threshold between model states showing 
enhanced and impaired selection remains (Supplementary results, 
Equation 24). Consequently, even with feed-forward inhibition of the 
response-representing populations, BHC-like changes always lead to 
impaired selection, and Huntington's disease-like changes can lead to 
enhanced selection, again providing evidence for our two hypotheses. 

3.7. Trade-off of NMDA sensitivity and connectivity loss is robust to 
nonlinear neuronal dynamics 

Such reduced models allow us a deep understanding of the micro- 
circuit's dynamics, but cannot account for the potential contributions 
of non-linearities in single neuron dynamics. Three non-linear pro- 
cesses are particularly relevant for our understanding of the striatal 
microcircuit in Huntington's disease and BHC conditions: the non- 
linear voltage-dependent activation and the long-lasting saturation 
of the NMDA receptors, and the reversal potential of GABAa receptors. 

As Huntington's disease is thought to lead to over-sensitive NMDA 
receptors (Fan and Raymond, 2007; Okamoto et al., 2009; Milnerwood 
et al., 2010), the first two processes are directly involved in how ex- 
citability of the medium spiny neuron changes in Huntington's dis- 
ease. Their voltage-dependence means that the NMDA receptors be- 
come increasingly available for activation as the neuron's membrane 
potential is increased towards its firing threshold (Jahr and Stevens, 
1990); this process feeds back on itself, whereby activation of some 
NMDA receptors depolarizes the neuron, which in turn makes further 
NMDA receptors available. In this way, NMDA receptors can non- 
linearly amplify the effects of small depolarizing inputs. The satu- 
ration of the NMDA receptors provides a ceiling on activity: as the 
neurotransmitter blocks the receptor for around 100 ms, so the re- 
ceptor is unavailable for further activation by incoming spikes from 
pre-synaptic sources (Lester et al., 1990). In this way, NMDA receptors 
can non-linearly reduce the effects of changes at high pre-synaptic 
firing rates. 

The reversal potential of GABAa synapses impacts the ability to 
perform selection, as the GABAergic connections between medium 
spiny neurons are not simply providing inhibition to each other 
(Plenz, 2003). The GABAa reversal potential is typically -60 mV and 
so is above the medium spiny neuron resting potential (typically -80 
to -90 mV) (Plenz, 2003; Bracci and Panzeri, 2006). Consequently, 
GABAergic input from other medium spiny neurons and fast spiking 
interneurons produces excitatory potentials in medium spiny neu- 
rons at rest and only produces inhibitory effects once the medium 
spiny neuron is sufficiently depolarized. 
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To examine the potential contribution of these non-linear phe- 
nomena, we constructed a new biophysical model of the striatal mi- 
crocircuit that incorporated a spiking medium-spiny neuron model 
with AM PA, NMDA, and GABAa receptors (Humphries et al., 2009a). 
The microcircuit model incorporated all three inhibitory connections 
examined in the reduced model: feedback inhibition between the 
medium spiny neuron populations, self-inhibition of the medium 
spiny neuron populations, and feed-forward inhibition from a fast 
spiking interneuron population (Fig. 4A). 

To isolate the contribution of single neuron dynamics from con- 
nectivity effects, we used a single model neuron per medium spiny 
neuron population, and modelled all other connection as sources of 
spike train inputs. We simulated the increased NMDA sensitivity of 
medium spiny neurons in Huntington's disease by increasing the 
NMDA receptor conductance in the medium spiny neuron models, 
and simulated medium spiny neuron loss (in Huntington's disease) 
and connectivity loss (in BHC) by reducing the number of inputs re- 
ceived from the global medium spiny neuron population. 

We found that this model reproduced the signal selection ability 
of the reduced model: a sufficiently large step in input to one popula- 
tion could cause that population's output to reach an upper threshold 
and the other population's output to reach a lower threshold. As the 
density of connections between the medium spiny neurons and from 
the fast spiking interneurons to the medium spiny neurons was de- 
rived from a detailed three-dimensional model of striatal anatomy 
(Humphries et al., 2010), and the synaptic conductances of these con- 
nections were taken directly from data (Koos et al., 2004; Gittis et 
al., 2010), this finding supports the above assertion that many, weak 
medium spiny neuron connections can produce signal selection. 

We found that our biophysical model qualitatively reproduced the 
results of the reduced model (Fig. 4B): increased NMDA sensitivity 
reduced the minimum input step needed for unambiguous response 
selection, whereas decreased medium spiny neuron connectivity gen- 
erally increased the minimum input step needed for unambiguous 
selection. Moreover, there existed two model states, one with en- 
hanced selection through sufficiently increased NMDA sensitivity to 
compensate for loss of medium spiny neuron connectivity, and one 
with impaired selection. Thus, our hypotheses for striatal competi- 
tion contributing to response selection and for increased NMDA sen- 
sitivity leading to enhanced response selection are not dependent on 
assumptions of linearity of the medium spiny neuron's dynamics. 



3.8. Selective vulnerability of indirect pathway striatal neuron 
populations still leads to trade-off 

A hallmark of the progression of Huntington's disease is that the 
D2-receptor expressing medium spiny neurons projecting to the ex- 
ternal pallidum are lost more rapidly than the Dl-receptor expressing 
medium spiny neurons projecting to the internal pallidum (Mitchell 
et al., 1999; Deng et al., 2004; Douaud et al., 2009). We thus also 
studied a reduced model incorporating this selective vulnerability. 
Each response-representing medium spiny neuron population was 
separated into its Dl and D2 receptor expressing sub-populations 
(Fig. 5A). Signal selection was read-out from the two Dl receptor 
sub-populations, as their target nuclei are thought to directly ex- 
press selection of motor programs (Mink, 1996; Redgrave et al., 1999; 
Hikosaka et al., 2000). 

Unlike the prior reduced model results (Fig. 3), analytical solutions 
are not available for this four-dimensional model and so numerical 
simulation was performed to find the minimum input step sufficient 
for unambiguous selection (Materials and methods). 

To mimic the selective vulnerability of D2 medium spiny neurons 
in Huntington's disease, we looked at the effect on selection of re- 
ducing only the inhibitory weights originating from the D2 medium 
spiny neuron populations, whilst increasing the cortical input weight 



to both types of medium spiny neuron. This was intentionally an ex- 
treme characterization of the selective vulnerability, with no loss of 
Dl medium spiny neurons at all: we sought to check whether or not 
loss of neurons from a sub-population would have the reverse effect 
of non-selectively losing neurons from the whole population. 

Fig. 5A shows that selectively losing only D2 medium spiny neu- 
rons increases the minimum input step required for selection, sim- 
ilar to the non-selective loss of medium spiny neurons. Moreover, 
Fig. 5A shows that correspondingly increased NMDA sensitivity in 
both Dl and D2 medium spiny neuron populations can still lead 
to the critical threshold between enhanced and impaired selection 
model states (white line in Fig. 5A). Neither the further addition of 
feedback self-inhibition for each sub-population (Fig. 5B) nor feed- 
forward inhibition from a fast spiking interneuron population to each 
sub-population (Fig. 5C) qualitatively changed these results. Thus, 
irrespective of the details of the microcircuit, if D2 medium spiny 
neurons are selectively vulnerable in Huntington's disease then their 
loss could still lead to impaired selection, which can be reversed by 
increased NMDA sensitivity. 

4. Discussion 

We tested patients with neurodegenerative disorders, including 
a rare case of BHC, against age-matched controls to dissociate the 
contributions of changes in the striatum to impaired and enhanced 
cognitive processes. Both the Huntington's disease patient group and 
the BHC patient had altered response selection on trials with irrele- 
vant sensory information, providing further evidence that changes in 
the striatum affect selection (Beste et al., 2008). The disorders caused 
opposite changes in response selection compared to controls: reduced 
errors and response time for Huntington's disease patients (replicat- 
ing our prior study, Beste et al., 2008); increased errors and response 
time for the BHC patient. Such deficits in cognitive functions in BHC 
are also shown in the neuropsychological test data (Stroop data) and 
are well in line with other results on cognitive functions in BHC (e.g. 
Gras et al., 2012). 

In the electrophysiological data, these behavioural differences 
were reflected in a more negative MMN and RON for Huntington's 
disease patients (replicating our prior study, Beste et al., 2008), and a 
reduced MMN and RON for the BHC patient. The P3a did not show 
modulation across the examined conditions, which suggests that 
frontostriatal processes do not affect involuntary attentional shifting 
mechanisms in this paradigm. 

Our results provide evidence that Huntington's disease and BHC 
have opposite effects on neural mechanisms for response selection: 
when medium spiny neuron dysfunction is combined with exagger- 
ated glutamatergic neural transmission (as in manifest Huntington's 
disease), the efficacy of response selection increases. When medium 
spiny neuron dysfunction alone is evident (as in BHC), performance 
declines. 

4.1. Dissociated contributions of NMDA sensitivity and medium spiny 
neuron dysfunction to selection 

We proposed that this dissociation between Huntington's disease 
and BHC effects on response selection supports two hypotheses. First, 
that computational processes involving the striatum directly influ- 
ence response selection in this task, and so changes in striatum are 
the cause of observed changes in selection. This is consistent with 
existing models and theories proposing that the striatum and the 
wider basal ganglia circuit underpin a general mechanism for selec- 
tion (Mink, 1996; Redgrave et al., 1999; Plenz, 2003). Second, that 
the enhanced selection seen in Huntington's disease patients is solely 
due to enhanced NMDA sensitivity compensating for medium spiny 
neuron dysfunction. 
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Our computational model showed that increasing NMDA sensitiv- 
ity indeed improved response selection compared to a healthy state 
model, whereas loss of medium spiny neuron local connections im- 
paired response selection. Crucially, the model also showed how these 
changes to response selection account for the performance changes 
of both disorders: loss of medium spiny neuron local connections 
alone, as in BHC, only impaired selection; yet many combinations of 
increased NMDA sensitivity with consequent loss of medium spiny 
neuron local connections, as in Huntington's disease, improved re- 
sponse selection. Our model thus predicts a critical threshold at which 
increased NMDA sensitivity can more than compensate for the loss 
of medium spiny neurons in Huntington's disease, resulting in the 
observed enhanced cognitive processes. 

We found that these results were robust to the details of the stri- 
atal microcircuit, to the selective vulnerability of the D2 medium 
spiny neuron population in Huntington's disease, and to the non- 
linear properties of NMDA and GABAa receptors. An open question 
is why the difference between groups was observed in the distrac- 
tor trials only. A key insight from our analytical results was that the 
minimum input step needed for selection increased in the order: "en- 
hanced" Huntington's disease-like state, the healthy state, and BHC- 
like state. Fig. SI shows how this ordering of input steps provides a 
plausible hypothesis for the error rate changes both between groups 
and between trial types. 

4.2. Enhanced or impaired selection is dependent on trade-off between 
sensitivity and cell loss 

Importantly, our models all suggest that enhancement of selection 
only occurred in a sub-set of all possible Huntington's disease-like 
conditions. Our analytical results show that whether Huntington's 
disease-like models show enhanced or impaired response selection 
depends on whether or not NMDA sensitivity increases beyond a criti- 
cal threshold for compensating the impairing effects of medium spiny 
neuron loss. Below this threshold, selection is worse than that for the 
healthy-state model. Our models thus predict that the normal im- 
pairment of Huntington's disease patients on a wide range decision- 
making tasks (Lawrence et al., 1998, 2000; Stout et al., 2001) is a 
consequence of NMDA receptor sensitivity not reaching this thresh- 
old. 

As a corollary of this result, our model predicts that our simple au- 
ditory task engages a region of the striatum in the "enhanced" state, 
where medium spiny neuron cell loss is not too advanced compared 
to increased NMDA sensitivity. According to the excito-toxicity hy- 
pothesis (Fan and Raymond, 2007) the two processes are causally 
linked: the excessive exposure to glutamate underlying the increased 
NMDA sensitivity results in the death of medium spiny neurons. We 
would thus expect that there would be an upper limit on the increase 
in NMDA sensitivity before its consequent loss of medium spiny neu- 
rons dominates. Indeed, the existence of a dorsal-ventral gradient of 
striatal cell loss in Huntington's disease (Douaud et al., 2009) sug- 
gests that a gradient of trade-offs between sensitivity increases and 
cell loss is plausible. The region of the striatum supporting the en- 
hanced performance is thus unlikely to be large: one possibility is 
that the auditoiy striatum plays a critical role in response selection 
in this task (Znamenskiy and Zador, 2013) and this region is in the 
"enhanced" state in Huntington's disease. 

It is conceivable that effects are not restricted to the auditory stria- 
tum, because NMDA neural transmission and medium spiny neuron 
density have until now not been shown to differ between striatal 
subregions. In this regard recent molecular genetic association stud- 
ies show that visual sensory memory is modulated by genetic poly- 
morphisms encoding for neurotrophic factors and receptors highly 
expressed in the striatum and of relevance to glutamatergic neu- 
ral transmission (e.g. Arning et al., 2014;Beste et al., 2012, 2011; 
Getzmann et al., 2013). However, we show that there could exist 



an "enhanced" selection state due to NMDA sensitivity compensating 
for cell-loss. But most studies show poor performance of Huntington's 
disease patients on response selection tasks. Our results suggest that 
in Huntington's disease most regions of the striatum are in the "im- 
paired" state, leading to the poor performance, but that the auditory 
task engages a region of striatum that is in the enhanced state — and 
that region is not usually engaged in cognitive tasks, otherwise we 
would have seen this enhancement of performance before. 

4.3. Potential limitations of the study 

However, it has to be acknowledged that Huntington's disease is 
also accompanied by changes at the cortical level (e.g. Tabrizi et al., 
2009), which was not taken into account in the model and which 
is clearly a limitation of the study. Obviously, these neocortical al- 
terations may further modulate the results obtained. As stated in 
the Introduction section, TITF1 has no known role in NMDA recep- 
tor function (Sussel et al., 1999). Whilst this may be regarded as a 
weakness in the a priori assumption of the presented study and the 
computational model (i.e., that there is no evidence falsifying an affec- 
tion of the NMDA receptor system in BHC), the experimental results 
presented for the Huntington's disease patients support the impor- 
tance of the NMDA receptor system for sensory memoiy and per- 
ceptual decision making processes in line with other data (e.g. Javitt 
et al., 1996; Kreitschmann-Andermahr et al., 2001; Umbricht et al., 
2002). The model simulations show that enhanced NMDA receptor 
involvement is critical for the effects to emerge. This, together with 
the fact that the BHC case and the Huntington's disease group reveal 
a completely dissociated pattern of results, replicating the pattern of 
a previous study (Beste et al., 2008), strongly suggests that NMDA 
receptor neurotransmission is likely to be unaltered in BHC. 

4.4. Further perspectives on decision-making 

In summary, we chose to focus here on understanding the putative 
striatal mechanisms underlying the observed differences in response 
selection between the control, Huntington's disease and BHC groups. 
This is a natural choice as the two disorders primarily affect the stria- 
tum, but it is likely that these experimental results will also provide 
fruitful ideas from other perspectives on decision-making and selec- 
tion. We have focussed on response selection here, where once a per- 
ceptual decision is made that decision is mapped to the appropriate 
response. Recent theories have suggested that the striatum and the 
wider basal ganglia circuit also play a role in the perceptual compo- 
nent of decision-making, particularly in the process of accumulating 
evidence for each hypothesized perceptual alternative (Bogacz and 
Gurney, 2007; Lepora and Gurney, 2012; Ding and Gold, 2013).Ding 
and Gold (2010, 2012, 2013) have provided intriguing evidence for 
striatal activity encoding both response and evidence accumulation 
components in separate neural populations, and that this activity 
causally affects the final decision. Consequently, it is possible that the 
changes to striatum occurring in Huntington's disease and BHC could 
also influence the perceptual decision (here, whether the tone was 
short or long) as well as response selection. However, standard mod- 
els of perceptual decision-making predict a speed-accuracy trade- 
off (Bogacz et al., 2006), whereby increased errors are the result of 
faster reaction times. Our data clearly show that both control and BHC 
groups increased both errors and reaction times on distractor trials: 
they thus present an interesting challenge to the models of perceptual 
decision-making. 
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